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Abstract 



We treat a tunneling electron coupled to acoustical phonons through a 
realistic electron phonon interaction: deformation potential and piezoelectric, 
in two or three-dimensional tunneling configurations. Making use of slowness 
of the phonon system compared to electron tunneling, and using a Green func- 
tion method for imaginary time, we are able to calculate the change in the 
transition probability due to the coupling to phonons. It is shown using stan- 
dard renormalization procedure that, contrary to the one-dimensional case, 
second order perturbation theory is sufficient in order to treat the deforma- 
tion potential coupling, which leads to a small correction to the transmission 
coefficient prefactor. In the case of piezoelectric coupling, which is found 
to be closely related to the piezoelectric polaron problem, vertex corrections 
need to be considered. Summing leading logarithmic terms, we show that 
the piezoelectric coupling leads to a significant change of the transmission 
coefficient. 

I. INTRODUCTION 

The effect of a bosonic environment on a tunneling particle has been extensively studied. 
The work of Caldeira and Leggett which originally studied the tunneling of a particle 
from a metastable state, focused attention on this problem. Today there are many separate 
fields devoted to related questions. As an example the case of coupling of a tunneling 
particle to a photonic environment is the basis of three such fields: that of dynamic image 
potential Coulomb blockade and dissipative macroscopic quantum tunneling (MQT) 
of Josephson phase || . The development of experimental techniques through which different 
effects now become within the scope of experimental verification || is a source of enhanced 
interest in such problems. 

The coupling of a tunneling electron to phonons has received much less attention p], 
though this problem is of fundamental interest, being closely related to the polaron problem 
0, as well as being of potentially great technological importance. In a previous paper || 
we have shown that in the one-dimensional case, i.e., a quantum wire or narrow constric- 
tion, the coupling to the zero-point fluctuations of the acoustic phonon field leads to an 
exponential enhancement of the tunneling probability. The question of what happens to 
this enhancement in the the two-dimensional case in which the electron's wave function is 
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restricted in one of the the directions perpendicular to the tunneling direction, and the three- 
dimensional case where the electron is free to move in both the perpendicular directions, is 
the subject of this paper. We address the problem of an electron tunneling coupled to acous- 
tic phonons via two realistic coupling mechanisms: deformation potential and piezoelectric. 
For both coupling mechanisms we calculated the change of the tunneling probability and 
the renormalization of the potential barrier height. In two and three dimensions the change 
of tunneling probability is smaller than in one dimensional case and for its calculation we 
developed a perturbation theory employing the Green function method. The problem of cal- 
culation of the phonon correction, however, appeared to be nontrivial for the piezoelectric 
coupling, where the summation of high order terms of the perturbation theory is necessary. 

Being interested in the effect of zero-point fluctuations of the phonon field on the tun- 
neling electron we did not consider the coupling to optical phonons. The virtual process of 
phonon emission and absorption is inversely proportional to the energy difference involved 
in that process. Since emission and absorption of optical phonons involves a large energy 
difference the interaction with optical phonons can be neglected. 

To calculate the transmission coefficient we make use of an approximation based on 
the slowness of the phonon field compared to the tunneling electron, which we refer to as 
the static approximation j5|, and obtain a stationary Schrodinger equation for the electron 
coupled to phonons. We then apply the parabolic (paraxial) approximation which is equiv- 
alent to considering the following physical situation (which is typical experimentally) (i) the 
barrier is so high that the tunneling can be considered semi-classically and (ii) the interac- 
tion energy between electrons and phonons is small compared to the height of the barrier. 
Treating the phonons via a Green function method the self-energy term is then calculated 
using second order perturbation theory. The resulting self-energy term is divergent at large 
transferred momenta, for both coupling mechanisms, and these divergences are eliminated 
by a renormalization procedure. For the deformation potential one has to renormalize the 
potential barrier height and the electron's effective mass, this finally leads to a correction 
to the transition probability. On the other hand for the piezoelectric coupling, after the 
renormalization of the potential barrier height the second order correction diverges loga- 
rithmically. The summation of main logarithmic terms resulted in a factor non-trivially 
depending on the electron-phonon coupling constant and can lead to a significantly change 
the transmission coefficient. 

The paper is composed in the following way. In Sec .[II] we formulate the problem, in- 
troduce the model, and present the static approximation which helps us facilitate the cal- 
culations. In Sec.|T| the Green function method through which we approach the problem 
is described, the parabolic approximation is introduced, self-energy diagrams are presented 
and it is shown how one uses them to derive the transition probability. In Sec.|y] the renor- 
malized self-energies for the two and three-dimensional tunneling configurations are given 
for the deformation potential coupling and the renormalization equations are presented. 
The resulting change in the transmission coefficient is calculated. In Sec.[V] the method for 
including vertex corrections for the piezoelectric case is given and the resulting change in 
the transmission coefficient is calculated. The calculation of some complicated integrals is 
given in appendices. 
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II. FORMULATION OF THE PROBLEM 



A. The model 



We treat the problem of an electron coupled to acoustic phonons while tunneling through 
a rectangular barrier, of width L. We will treat two cases that of a two and three-dimensional 
tunneling configurations. In the three-dimensional, case the electron is free to move in 
the perpendicular (to the tunneling) directions, whereas in the two-dimensional case the 
electron's movement is restricted in one of the perpendicular directions i.e a two dimensional 
electron gas. In both cases phonons move freely through the bulk and thus are treated as 
three-dimensional. Being interested in the effect of the zero-point fluctuations of the phonon 
field on the electron tunneling probability we consider the system at zero temperature. 



In the general case of a tunneling particle coupled to the environment two typical time 
scales arise: the traverse time for the tunneling particle and the typical time in which 
the environment reacts to the tunneling of the particle. When the tunneling particle is 
an electron and the environment acoustical phonons, it can be shown that typically the 
traverse time for the electron is much smaller than the environmental reaction time scale. 
The small parameter expressing the small ratio between these time-scales is: s/v , where s 
is the phonon sound velocity, and v is the under the barrier electron semi-classical velocity. 
|| For a potential barrier of around 10 meV in bulk GaAs, this ratio is of the order of 
0.01. This small ratio leads us to what we refer to as the static approximation. In this 
approximation the phonon dynamics is neglected and the potential created by the barrier 
has a phonon addition that may have different time independent shapes with the probability 
equal the wave function squared of the phonon field. [For a more detailed description of the 
static approximation one is referred to the previous paper ||.] As a result the electron wave 
function underneath the barrier satisfies the following stationary Schrodinger equation for 
the electron coupled to phonons 



where x is the electron coordinate, m is the electrons effective mass, V is the potential 
barrier height, u q = u_ q * are the Fourier coefficients of the lattice displacement vector 
«(x). u q is a random variable whose probability is controlled by the ground state phonon 
wave function. For the two-dimensional case the Laplacian is the two-dimensional Laplacian 
operator, /(x, w q ) is the interaction term between the electron and phonons defined in the 
following way. For the three-dimensional case, in which the electron wave function is not 
restricted in the perpendicular to tunneling direction which is taken as x, the deformation 
potential coupling interaction term is given by 



B. Static approximation 




(2.1) 




(2.2) 
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where A is the deformation potential constant and V vo i is the normalization volume. For the 
piezoelectric coupling the interaction term in the Hamiltonian is given by 



/p(x,u q ) = —. ==^w q exp(zqx), (2.3) 

V V vo l q 



where H = 47rey (Pijk 2 )n/^, e is the electron charge, e is the dielectric constant, is the 
piezoelectric modules and denotes the average over angles. Though the piezoelectric 

interaction is anisotropic the isotropic Hamiltonian approximation is commonly used |5j. 
Subscripts ti and p shall be used throughout this paper, unless otherwise stated to define 
deformation potential and piezoelectric related terms respectively. 

For two dimensional electron gas were the electron wave function is taken as restricted 
in the z direction the interaction terms are given by 

iA 

/ d (x, u q ) = —= |q| Mq M qz exp [i(q x x + q y y)\, (2.4) 

V "vol q 

for the deformation potential coupling and 

/p(x,u q ) = -7= 51 M q M 9z exp [i(gr x x + g„y)]. (2.5) 

V ^uoi q 

for the piezoelectric coupling. Here M qz is the matrix element of the phonon exponent, 
exp (iq z z) between the wave functions describing the electron quantization in the z direction. 

We neglect screening of the electron-phonon interaction, since underneath the potential 
barrier there are no free electrons and screening by remote electrons is small. 

Although the static approximation resembles the electron tunneling in the presence of 
static disorder there is a very important difference between the two problems. Whereas a 
static random field can absorb any momentum from the electron the momentum, absorption 
in the static approximation is related to phonon emission. Since under the static approxi- 
mation it is assumed that the phonon configuration does not change the tunneling electron 
momentum is conserved. 



III. GREEN FUNCTION METHOD 

The commonly used approach for addressing the problem of a tunneling particle coupled 
to the environment is based on the path-integral approach. The path-integral approach to 
the problem is extremely convenient in the one-dimensional case where the contribution to 
the transition probability is essentially from one trajectory. In two and three dimensional 
cases many different trajectories contribute to the transition probability and we prefer to 
use Green function method. The application of this method is facilitated by the fact that a 
reasonable potential barrier (e.g., 10 meV) it is much larger than electron-phonon interaction, 
so for the calculation of electron Green function we can use perturbation theory. In some 
points the Green function method that we use differs from commonly used formalism and 
we present it's derivation. 
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A. Parabolic approximation 



We consider the case when the coupling between the electron and the acoustic phonons 
is small compared to the effective potential barrier height V — E, which is considered as 
large. More than this, we assume that h/^j2m(V — E) is the smallest length scale in the 
problem. It is much smaller than the width of the barrier L and in two-dimensional case 
it is much smaller than the size of the quantization in z direction, d. The first of these 
conditions enables us to use the parabolic (paraxial) approximation |10| , and is the basis for 
our perturbative method. 

Under the parabolic approximation we search for a solution to the Schrodinger equa- 



tion (O) of the form: ^(x, u q ) = exp(—Kx)ifj(x } Uq) where x is the tunneling direc- 



tion, the exponential factor characterizes tunneling without interaction with phonons, 



k = J2m(V — E)/h, and ip(x,u q ) is a slowly varying envelope function. Slow variation 
of i/j(x,u q ) follows from the weakness of electron-phonon interaction and leads to the pos- 
sibility of neglect of its second derivative with respect to x, (d 2 ip /dx 2 ) <C n(dip/dx). Then 
Eq. ( |2.1| ) becomes 

| 011) TFl 

where Vj_ designates the gradient with respect to the perpendicular to tunneling direction 
x (in the two-dimensional case Vj_ = d/dy). This Schrodinger equation is just equivalent 
to a non-stationary Schrodinger equation where x/k plays the role of imaginary time and 
there are two spatial dimensions instead of the original three. 

The transmission across the barrier is characterized by the solution to Eq. (|3.1| ) with the 
boundary ("initial") condition ip(x = 0, rj_) = 1. It can be written as 

ip(x, rjwq) = J g(x, r±; 0, r'Jwq^O, r' ± )dr' ± . (3.2) 

Then the transmission coefficient is given by (multiplied by exp(— 2kV) which is due to the 
parabolic approximation) 

|V(L,rJ Mq )| 2 ) = J (G(L, r_L] 0, r^|w q )^*(L, rj_; 0, r" |« q )) -0(0, r^)-0(O, r'[)dr' ± dr'[ , (3.3) 

where the angular brackets mean the averaging over the phonon configuration. Due to 
translational invariance the average of the Green function product depend on the difference 
r'j_ — r" so that the rhs does not depend of rj_. 

In the case of tunneling in the presence of a random field the average of the Green function 
product mixes both Green functions, i.e., it is not equal to the product of averages. The 
situation is much simpler in the case of phonons under the static approximation. The neglect 
of phonon emission means that ip(L, r±\uq) in the lhs of Eq. ( |3.3| ) corresponds to the phonon 
ground state. In other words, for the calculation of the transmission coefficient we have to 
keep only the diagonal with respect to phonons part of the Green function. This diagonal 
part is obtained by the averaging of Q over the phonon configuration, G(x, rj_ — r'j_) = 
(£(x,r ± ;0,r'Jw q )). 
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It is convenient to calculate G(x, r± — r' ± ) in the Fourier representation. Actually we 
perform Fourier transformation for the perpendicular direction and a Laplace transformation 
for the x direction. After these transformations Eq.( |3.1| ) becomes 

1 

-fclV'kxb) + k5u ± ,o + Kpip k± (p) = EE F q±(<7z)^k ± -q x (p - iq x ), (3.4) 

qx qx 

where ipk ± (p) are the Fourier and Laplace transforms of u q ) with respect to the perpen- 
dicular and x directions respectively and F q± (g x .) = — (m/h 2 )/^^). Throughout the paper 
we will use the notation p for the Laplace transformation parameter (imaginary x-component 
of the electron wave vector), for the electron wave vector in the plane perpendicular to 
the tunneling direction and q for the phonon wave vector. 
The transmission coefficient now is \t\ 2 e~ 2K-L where 

t = — / dp e pL G k±=0 (p) . (3.5) 

Calculating the transmission coefficient thus amounts to calculating the averaged Green 
function. 

B. Green function calculation 

The Green function unperturbed by the phonon potential is G^ Op) = (up + k\/2)~ l so 

(3.6) 

where S k± (p) is the self energy. 

We calculate the self-energy with the help of a diagrammatic technique that differs from 
the regular one in random field by two points. The first is that the energy variable is 
replaced by — up. The other point, that in the static approximation the two-particle Green 
function, Eq.( [3.3|) , equals the product of two one-particle Green functions, has already been 
mentioned. It is also possible to demonstrate this last point diagrammaticly, showing that 
in the static approximation, due to momentum conservation, diagrams with random field 
lines connecting two Green functions contain the factor l/V, where V is the volume of the 
system. 

Due to the weak coupling between electrons and phonons, in most of the cases that we 
consider a good approximation of the self-energy is given by the first order diagram presented 
in Fig.||, i.e., 

£k» = E (I^X-qxfe)lX >-^)> ( 3 - 7 ) 

The averaging over phonons is carried out with the help of the phonon wave function 
squared, C exp(— X) q p|q|s|w q | 2 /^), where p is the the crystal density s is the sound velocity, 
and C is the normalization constant. Since F <i± (q x ) is proportional to w q the averaging gives 
(F q± F q / ) = a|q , | ±1 5 qX) _ q ' . The +1 corresponding to the deformation potential coupling 



GkAp)= KP + k 2 ± /2-E k± (p) 
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and the —1 sign corresponding to the piezoelectric coupling. For the deformation potential 
coupling a is a dimension-less parameter given by 

a d = • (3.8) 

n /m z 

For the piezoelectric coupling mechanism a has units of one over length squared and is 
defined as 



h/ps _ 2 

h 4 /m 2 

For the calculation of the transmission coefficient (|3.5| ) we need the self-energy only for 

e *x=o(p) = (gp £, Jfc / <*V <p _ ^ + i q±2 (Vol 3 + <z* 2 f\ (3-10) 

the (d, +1) and (p, — 1) subscripts and signs correspond to the deformation potential and 
piezoelectric coupling respectively. 

The self energy for the two-dimensional tunneling configuration is given by 

E 2 (p) = W^T dq z \M q f r dq x r dq y - L-_ + g 2 + 2 ) ± \ (3 . n) 

where the subscript 2 will be used to define the two-dimensional case throughout the paper. 



IV. DEFORMATION POTENTIAL COUPLING 

The self-energy integrals for the three-dimensional tunneling configuration (|3.10|) have 
an ultra-violet divergence for both coupling mechanisms. Actually these divergences are 
cut off by the inverse lattice constant a, but it goes against physical intuition when for a 
phenomena involving momenta much smaller than the inverse lattice constant such cut-off is 
necessary. The divergences can be eliminated with the help of the standard renormalization 
procedure (see, e.g., Ref. [0). For the deformation potential coupling the large constant 
term S(0) (which has a form ad(co/a 2 ) where Co is a numerical constant) can be viewed as 
a renormalization of the potential barrier height, V and, hence, the renormalization of the 
propagation constant k, 

K 2 ren = 2m(V -E) + S d (0). (4.1) 

The abbreviation ren is used to denote renormalized quantities, also it is assumed that k is 
K ren (unless otherwise stated), through the rest of the paper. Actually the renormalization 
Q4.1D has to be made before one passes to the parabolic approximation. 

Canceling the non-p dependent diverging part of for the case of deformation poten- 
tial coupling does not cancel all the divergences. One needs further to cancel a term S^(0)p, 
(Srf(0) is logarithmically divergent, it can be estimated as a^KCi ln(/ta) where C\ is a numer- 
ical constant) which can be viewed as a further renormalization due to a renormalization of 
the electron's effective mass, 
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S d (0) 



77? 



(4.2) 



Subtracting the two diverging terms from the self-energy integral one obtains 

K 2 \ 3 



= s d (p) - S d (0) - S d (0)p 



-«d- 



67T 2 



ln2A 



(4.3) 



where A = p/k. 

For the two-dimensional tunneling in the deformation potential coupling the self-energy 
also has an ultra-violet divergence. After subtracting the divergent term from the self-energy 
integral, we obtain a finite expression, for the self-energy 



E2T(p) 



£ rf , 2 (p) - s d,2(0) = -a d 



3/2 



2tt 2 



-P 



dq z \M q 



(4.4) 



For details of this calculation one is referred to the Appendix [A|. 

In the tunneling problem, contrary to the propagation one, in the case of weak coupling 
the Green function can be expanded in the self-energy. Then after the substitution of Eqs. 
( |4.3| ) and ( O ) into Eq. ( |3.5| ) the inverse Laplace transformation for the three-dimensional 
case one gets that — 1 — (ad/Q^ 2 ){nL)~ 2 . Thus the transition probability for the electron 
coupled via the deformation potential coupling mechanism is given by 



rp rpO 



3tt 2 ( K Ly 



(4.5) 



where T° is the transition probability for the tunneling problem without phonons. 
In the two-dimensional tunneling case the transition probability is given by 



rpO 

1 2 



1 - Ol d 



3/2-1 

TT 2 K 



dq z \M g 



(4.6) 



V. PIEZOELECTRIC COUPLING 

In the two-dimensional case the integral for the piezoelectric self-energy is finite. The 
finite independent of p part of the integral, S p (0) ~ (ap/dK)\M \ 2 \n(dK), can be viewed as 
a finite renormalization of k, where M is the matrix element of the phonon exponent, for 
q z = 0, and d is the length of the potential barrier in the z direction. 1/d is the scale that 
cuts off the q z integration, and it is assumed that in the physical situation we are studying, 
1 dn. Renormalizing E we have 

^» = Sp, 2 (p) " Sp, 2 (0) = 7^ \ |Mo| 2 In (pd) (5.1) 

In the three-dimensional case the renormalization correction diverges in the parabolic 
approximation and the renormalized self-energy is 



S 



E«» = E p (p) - E„(0) = ^ A (In A + ln2 - 1) . (5.2) 

Details of the calculation of the integrals can be found Appendix [A[ 

The terms A In A in Eq.( |5.2|) and pln(pd) in Eq.( |5.1D become large compared to up as p 
goes to zero. Thus the perturbation theory in this case is not sufficient. 

A similar problem was encountered by Edel'shtein when he studied the piezoelectric 
polaron ||12| . He considered high order vertex corrections and showed that in each order of 
the perturbation theory the highest power of the logarithm comes from diagrams with non- 
crossing phonon lines. Neglecting all other diagrams he succeeded in obtaining an equation 
for a vertex correction. Because the origin of the logarithmic divergence in Edel'shtein's 
case and ours is the same we can use the same approximation. In the following calculation 



we follow closely the derivation presented in Ref. [|12[ . 

In the following discussion and calculation we will treat the three dimensional case the 
treatment of the two dimensional case is very similar therefore only the results shall be 
presented for it. 

We start by considering vertex corrections. The first vertex correction, T\, presented in 
Fig- i 

r x (pi, ku_,p 2 , k 2 ±) = a p J GSx-aiCPi - z ^) G! k 2 ) x -q ± (P2 - %q*). (5.3) 

The integral is calculated in Appendix |B| giving 

where rj = up + k\j2 (p is either pi or p 2 , since under logarithmic accuracy the difference 
between them can be ignored, for the same reason k± is either ki± or &2±)- For k± = the 
same result can be obtained from the self-energy ( |5.2p with the help of the Ward identity. 
For the two dimensional case the first vertex correction is 

a P \ I i/ |2 • 



^-i^ji^n^j' (5 - 5) 

where fj = up + ky/2. The calculation of this integral can also be found in Appendix [B|. 

The second order diagrams for the vertex corrections are presented in Fig.[| The diagram 
(3, a) is reduced to the renormalization of the Green function. It is possible to show that 
the leading order contribution to the diagram (3,d) is of the order of (a p /27r 2 K 2 ) 2 ln(^//t 2 ) 
containing lower power of the logarithm than the power of a p . The leading order contribution 
to the other two diagram is of the order of (a p /2% 2 K 2 ) 2 ln 2 (r//fi: 2 ), i.e., the power of the 
logarithm equals the power of the coupling constant. The same relation holds true for 
more complicated diagrams. That is, diagrams containing intersections of phonon lines 
contain lower power of the logarithm than the power of the coupling constant and diagrams 
without such an intersection get another power of the logarithm in every next order of the 
perturbation theory. 

The last statement can be proven in the following way. The increase of the order of 
the perturbation theory by unity in a vertex diagram without intersection of phonon lines 
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corresponds to the replacement of one of the bare vertices in this diagram by IV Similar to 
the calculation in Appendix [B|, the values of the momentum in the contour of this new Ti 
that give the main contribution are much larger than the momentum in external contours 
but much smaller than k. So this Ti again gives a logarithmic contribution. 

Another result of such a calculation is that in diagrams without phonon line crossings 
the variables p and k± enter only in the combination 1] = up + k\/2. 

Thus when t]/k 2 is small enough, e.g. the effective potential barrier is high enough, 
one can neglect diagrams containing phonon intersection lines. The summation of all other 
diagrams leads to the equation for the vertex T presented in Fig.[|. 

Since the Green function and the vertex correction are functions of rj, the Ward identity 
is expressed in the following form: 



Integrating Eq. ( |5.6|) and making use of the fact T is a slow function we get 



lation of two vertices and resulting in a linear equation for the vertex function. The main 
contribution to the integral term in this equation comes from the region where the mo- 
mentum in the contour in Fig.|] is much larger than the external momenta. Similar to the 
calculation of I\, the difference between the external momenta is not important and they 
have to be kept only to limit the value of otherwise logarithmically divergent integral. So 
the equation can be written in the following form 



In this integral t ~ 1 and the terms q 2 /2 and r] are small compared to nq. So these terms 
can be neglected in the argument of T (it is a slow function) and t 2 can be neglected in the 
coefficient of q 2 /2 in the denominator (the exact value of this term is not important). Now it 
is convenient to consider T(r]) as a function of the large logarithmic variables x = \n(i]/ 'k 2 ). 
Then keeping just the large logarithm x' = ln(g//t) in the argument of T in the integrand 
we can easily integrate with respect to t. As a result, with the logarithmic accuracy, 





(5.8) 




(5.9) 



The solution of this equation and using Eq.( |5.9|) gives 




(5.10) 



In the two dimensional case the Green function is given by: 




(5.11) 
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Inserting these results in the equation for the transition amplitude Q3.5|) and performing 
the inverse Laplace transform [13| we obtain for the three-dimensional tunneling configura- 
tion 



T = T° (kL)~^ . 
The two-dimensional case gives the following transition probability 

felMol 2 



rp rpij 



I d\ **** 



L) \L / 



(5.12) 



(5.13) 



VI. SUMMARY AND CONCLUSIONS 



In this paper we have presented a detailed study of the effect of coupling of an electron 
tunneling across a rectangular barrier, to the zero-point fluctuations of the acoustic phonon 
field. Two coupling mechanisms, the piezoelectric and deformation potential, were studied 
for two and three dimensional tunneling configurations. We considered a small coupling and 
used two main approximations: the static approximation and the parabolic approximation. 

Contrary to the one-dimensional case 0, the main effect of the deformation potential 
coupling is the barrier and the electron mass renormalization. On the top of this the sec- 
ond order perturbation theory is sufficient. As a result only a prefactor in the tunneling 
probability is changes. 

For the piezoelectric coupling mechanism, the situation is quite different. The pertur- 
bation theory is not sufficient and an accurate calculation of the vertex is necessary. As a 
result of the coupling, in both two-dimensional and three-dimensional cases the tunneling 
probability acquires a prefactor proportional to a power of the barrier width. The small 
coupling constants enters in the power. This result resembles piezoelectric polaron. When 
the barrier is wide enough the piezoelectric coupling can lead to a substantial reduction of 
the tunneling probability. 

It is interesting also to compare the sign of the coupling effect in one- and higher di- 
mensions. In one-dimensional case the coupling led to an enhancement of the tunneling 
probability. The physics behind this is that an electron chooses for tunneling the most fa- 
vorable phonon configuration. || In high-dimensional case, due to ultra-violet divergences, 
such favorable phonon configurations led to overall renormalization of the barrier which 
cannot be separated from the bare barrier. The measurable effect depending on the length 
of the barrier is a reduction of the tunneling probability. 

In both cases the effects mentioned are enhanced when tunneling is near the top of the 
barrier, though as the effective potential barrier height, V — E, is lowered two effects become 
important: dynamical corrections to our static approximation and an effect which can be 
described as channeling, the change of the electron wave function due to the interaction with 
phonons. These are subjects for future study. 
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APPENDIX A: CALCULATION OF THE SELF-ENERGIES 



In this appendix we present the main stages in the calculation of the self-energy integrals 
appearing in the text, except the self-energy integral for the three-dimensional deformation 
potential coupling which is straight forwardly evaluated. We start with the calculation of 
the self-energy integral, for the deformation potential coupling in the two-dimensional case, 
after the subtraction of the divergent term. 



fy , poo POO POO 

^ en (p) = -^ / dq z \M g f dqj dq y 

Z7T J —oo J —oo J —oo 



a d f°° j i If ,2 f°° j f°° j ( 1 1 

1^2 ~ I 1„2 



(27r) J-oo z J-oo J-oo \Kp - iKq x + ~q% + f^/ 

x + <ly + ^ ■ ( Al ) 

Noticing that for the q x integration poles are located in the lower half of the complex q x 
plane we shift the q x integration contour to the upper half plane. T he shifted contour goes 
from ioo along the imaginary axis, rounds the point x = i\Jq% + q z in the positive direction 
and comes back to ioo. Thus we obtain 



S 2 e » = j^Jq z \M q f j™Jq y 



V^i dx W+^ + h 2 y q2y ql ' (A2) 



Where the upper cut-off, k, for the x integration is due to the parabolic approximation. The 
integration with respect to q z is cut off by the matrix element at q z ~ 1/d where d is the 
size of the quantization in the z direction. We assume that dn ^> 1 and hence the main 
contribution to the integral ( |A2| ) comes from the region q y ^> q z . Then we can neglect the 
q 2 z appearing in the square roots, and 



a d np 



W-Tf / dq z \M qz n dq y dx v (A3) 

A-k a J-oo J-oo J\q y \ [up + kx + ^){kx + 

Since we are interested in the limit p ~ 1/L which is much smaller than all other scales, we 
neglect up denominator, scaling out k by changing variables x to nx and q y to ny we arrive 
at 

^(P) = -S? f dq z \M g f \°°dy f'dx /f'zt • ( A4 ) 



Changing variables to: x = ty and changing the order of integration with respect to t and 
y we get 

s "» = - a s l d ^ r dt ^r jrrw • <A5) 



Were this integral is already easily evaluated resulting in Eq. 

We now evaluate the self energy term for the piezoelectric coupling in the two-dimensional 

case, 
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S^P) = 7^1 / dq z \M q f dqj dq y 

(z7Tj J— oo <J — oo J ~oo \ 



np - inq x + \q 2 -wq x + \q 2 y / 



x . (A6) 

7* + 9y + 9z 



Using the same contour integration we derive 

a p Kp r°° 2 f°° , t°° dx 



V /^|TS («p + kx + \qf){KX + 2<?2) 

x 1 . (A7) 

/x 2 -q 2 y -q 2 z 

q y in the integrand is smaller than x and as we will see iCk. Due to this we neglect g 2 
compared to kx and we obtain 

ST(P) = /°° AfcW fX /"^ , „ ^ n T^Z ■ (A8) 



x/i+S ^x 1 -q 2 y -q 2 z (P + ' 



Changing variables to: x = ty g 2 + g 2 ; = yq z , we get 

Er(P) = r iqAM j r * r * — * . 

27r k j-oo Jo + 1 -/i tvt - 1 p + *vy + l<?z 

When p <C 1/d the integral with respect to q z contains a large logarithm. In the logarithmic 
approximation 

arc?) = -^iM i 2 r- r^- r ^ , ( ai o) 



which results in Eq. (|5.1| ). 

Finally we calculate the self-energy for piezoelectric coupling in the three-dimensional 

case, 

Tr n (p) = [°° dq x ( d 2 q ± ( 1 —r^ I ; 1 (All) 

pW (2tt) 3 J-oo qx J HX \K V -iKq x + \ql -mq x + \q\) JqJ+ql 



Performing the angle integration and scaling out K and changing variables to: q]_ = 

2ttpz; q x = px, we get 

rv roc roc "I "I 

Tr en (p) = ^ / dz / dx- ; (A12) 

PK > (27T) 2 JO Loo (1-IX + Z)(-IX + Z) ^ X 2 + 2z/\ 

where A = p/ k. Noticing that for the x integration poles are located in the lower half of the 
complex x plane we shift the contour as we did before. Now the shifted contour rounds the 



point x = i\2z/X and the result becomes 
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■yren 



(P) 



a, 



2vr 2 



dz 



dy 



1 



a p X r°° 



2vr 2 Jo 



dt 



1 + ^f + z + y) ^f + z + y) J 2y 
du 



yj(2 + t)t J ° (l + U + Y + «*) (l + T + *) 



In the integral with respect to u we neglect terms of the order of A so that 



«pA r 00 



c/i 



In - - 1 + 21n(l +t) 
A 



a p A 
2^2 



In- - 1 + 2(1 -In 2) 
A 



+ y 2 



(A13) 



(A14) 



APPENDIX B: LEADING ORDER CALCULATION OF FIRST VERTEX 

CORRECTIONS 



Under logarithmic accuracy the main contribution to the integral ( |5.3| ) comes from the 
region k ^> q ^> p, h±. For this reason difference between pi and P2 as well as between ku_ 
and k 2 j_ is not important. So the integral (|573|) can be written as 



p / \ f°° qdq r dOsmO 

Jo (2tt)^o (-i qKC os9 + £ sin 2 9 + rA 



where rj = up + k 2 L /2. The term g 2 /2 cuts off the integration with respect to q from above 
and the term 77 does the same from below. In the region between these limits the integral is 
logarithmic and its value is accumulated from the whole interval of the integration. So that 
the exact value of these limits is not important and we can simplify the integration by the 
replacement of sin 2 9 with unity. As a result 

r f \ 2a P [°° ^dq Up . f V \ / n0 N 

TM = Wf^L n 2 7 (V 7 ^ 2= "w ln U*J- (B2) 

Where a lograthmic approximation has been used in calculating the last integral. 

Considering the first vertex correction in the two dimensional case again under loga- 
rithmic accuracy the difference between p 1 and p 2 as well as between k ly and k 2y is not 
important. Thus the first vertex correction is given by 

r i(^) =7?*r *z.mj 2 r r r= ^ — ^ — ^ > m 



(27T) 



where rj = np + k 2 /2. Using the contour integration described in Appendix (0) one gets an 
equation which is very similar to Eq. ( |A7|j 

ttp f°° , , , , ,2 Z" 00 . Z" 00 ^ 
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Treating this integral in the same way Eq. (|A7|) has been evaluated results in Eq. G5.5|). 
The main difference between this I\ and that for three dimensional case is that now the 
logarithmic integral is cut off from above by the upper limit of q z which is of the order of 
1/d < k. 
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FIGURES 



FIG. 1. The self energy term. The solid line represents the bare electron propagator whereas 
the dashed line represents interaction with phonons. 



FIG. 2. The first correction to the vertex T. 

FIG. 3. Second order vertex corrections 

FIG. 4. The equation for the renormalized vertex T resulting from the summation of leading 
logarithmic corrections. Double lines represent the full electron Green function. 
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